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Abstract 

It was recently pointed out that the decays D*^^ and D*^^ 

can be used for an extraction of \Vub\- The theory of these decays is poorly 
understood. It was shown that in a world of almost degenerate b and c-quarks 
the decay would be computable. The severe difficulties that are encountered in 
the realistic calculation stem primarily from the very hard photon produced 
in the two body decay. We point out that in the decays B^ — > D*^e^e~ 
and B~^ D*^e^e^ the photon vertex is soft when the charmed meson is 
nearly at rest (in the B~^ rest frame). This allows us to compute with some 
confidence the decay rate in a restricted but interesting kinematic regime. 
Given enough data the extraction of V^b with reasonably small uncertainties 
could proceed through an analysis of these exclusive decays much as is done 
in the determination of Vcb- 
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I. INTRODUCTION 



The determination of the Kobayashi-Maskawa (KM) matrix element Vub is one of the 
central themes of particle physics today. The verification of the KM theory of CP violation 
in the six quark model requires an independent measurement of the magnitude and phases 
of the KM matrix elements. The physics community has invested heavily in 5-factories, 
facilities devoted to the measurement of these magnitudes and phases. However, the extrac- 
tion of Vub from data requires calculation of hadronic transition rates that are marred by 
the usual host of difficulties introduced by the strong interactions. 

For example, the first determination of V^f, has been made in the inclusive semileptonic 
decay of i?-mesons. The inclusive decay rate, appropriately smeared over kinematic vari- 
ables, can be reliably estimated Since ^ iKifeP, the rate is dominated by the 
decay into charmed mesons. To isolate the contribution from Vub measurements are made 
of the decay rate for electron or muon energies that are large enough to exclude contam- 
ination from the decay into charm. However, in this very restricted kinematic range the 
theoretical prediction is highly unreliable, and it is known to be dominated by incalculable 
non-perturbative effects 0. 

A novel procedure offers more hope. If one can reconstruct the neutrino in the semilep- 
tonic decay of B mesons only small invariant hadronic masses are required to suppress or 
eliminate the contribution from decays into charmed particles 0. The invariant mass spec- 
tra are under better theoretical control because the allowed final states are more numerous, 
so quark-hadron duality is more likely to be valid than in the restricted region of very high 
electron or muon energies. 

In this paper we investigate the extraction of Vub from an exclusive decay mode of charged 
i?-mesons. Exclusive semileptonic decay modes, such a.s B ^ Triu or B ^ piu, have proved 
difficult to study theoretically. The rates are described by non-perturbative form factors 
which are difficult to estimate. The exclusive decay mode B'^ — > D*~^e~^e~ depends on non- 
perturbative dynamics too. But as we will see, under reasonable theoretical assumptions, 
the non-perturbative dynamics can be extracted from other measurements. 

A similar idea was recently proposed in the description of the decays B^ D*'^'-/ and 
B^ — s> D*^--f 0. There the transition was assumed to be dominated by the lowest poles {B* 
and Dg or D) and was given in terms of the amplitude for BB mixing and the probability 
rate for D* D'j. The validity of the symmetry arguments used there relied heavily on a 
theoretical limit, that rrib — rric ^ Aqcd, which is not likely to be a good approximation to 
reality. This assumption was important in relating the hard photon emission in B^ —>■ D*~^'y 
to the soft emission in D* —>■ D'j. We observe in this paper that for the decay into a lepton 
pair, rather than a photon, the extended phase space includes a region where the D* meson is 
nearly at rest in the B rest-frame, in which the photon is soft. Consequently the unphysical 
assumption rrib — rric ^ Aqcd can be dropped. 

We make several assumptions and approximations in our calculation of the rates for 
B^ — s> Dl^e^e~ and B^ D*^e^e^. If the assumptions prove to hold and the approxi- 
mations are reliable a determination of Vub through these exclusive modes could proceed in 
much the same fashion as that of Vcb through the exclusive semileptonic mode B D*£i'. 
Theory reliably predicts the rate in the kinematic region in which the D* is produced at rest 
in the B rest-frame. Experimentally the rate is measured away from this point and then 
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extrapolated. The procedure would be entirely analogous for the much rarer, but cleaner, 
decay we study, D*^e^e~ . 

II. BASIC ASSUMPTIONS AND SYMMETRIES 

Before we present our calculation we discuss in detail the nature of the assumptions and 
approximations that we use. To lowest order in electromagnetic and weak interactions, the 
amplitude for the process D*~^'y* is given in terms of 

{D*+\ J d'x e">-^ T{jitJx)O{0))\B+). (2.1) 

Here O = Jry'^ {l—'y^)uc'~f'^ {l—'~f5)d is a four quark operator responsible for the weak transition, 
j^jn is the electromagnetic current and T stands for time ordering of these operators. For 
each time ordering one may insert a complete set of states. Our first assumption is that the 
lowest mass states dominate, giving 

Jd<!> jd^x e*'?-^(^(x°) (D+l 0(0) 

+ e{-x'){D*^\OmB*^){B*^\3U^)\B^)\ (2.2) 

where / (i$ is the one body phase space. We know of no reliable way to test this assumption. 
Much of heavy quark chiral perturbation theory p relies on analogous assumptions. The 
effect of excited heavy mesons in loops in chiral perturbation theory has been analyzed and 
the result seems to be sensitive to the precise values of masses and couplings 0. This is an 
area where more research is needed before firm conclusions can be drawn. 

Given this assumption, the transition amplitude is calculated in terms of the matrix 
elements for the weak operator 

(D+|O(0)|B+) and (D*+|C(0)|5*+) (2.3) 

and the matrix elements of the electromagnetic current, 

{D*^\3U^)\D^) and (5*+br^(0)|i?+). (2.4) 

The matrix element of a vector current between a pseudoscalar and a vector meson can be 
written in terms of a single real form factor, which can be defined as follows 

(P', e|jr^(0)|p) = igxitV^^'PuP'^C- (2.5) 

Here X = B or D identifies the one particle states, e is the polarization of the vector 
meson, and the time component of the four vectors p and p' is the energy of the state, 
Po = \l M"^ + |pp. The Lorentz invariant form factor Qx is a function only of the invariant 
momentum transfer t = {p — p'Y ■ Heavy quark flavor symmetry |^ can be used to relate 
and Qjj when the argument is t ~ 0. Moreover, Qxi^) = /^X' where fix is the magnetic 
moment defined in Ref. |^ determines the radiative decay rate of the D* meson 
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The heavy quark contribution to the magnetic moment of the meson is suppressed in the 
large mass hmit, so fix = f^Qq, where Qg is the charge of the hght quark or anti-quark in 
the heavy meson. Therefore ~ \iJ'b+- 

The functional form of the form factor g{t) is not known. Clearly the dimensionful scale 
that determines its behavior is on the order of the mass of the lightest vector meson, the p- 
meson. Now, in the rest-frame of the decaying meson, p = 0, and t = (pq — M)^ — |p'p = 
{M' — My —2M {pq — M') . So at p' = the momentum transfer variable is small and positive, 
t ^ (50 MeV)^ for B mesons and t ^ (150 MeV)^ for D mesons. As |p'| increases t goes 
through zero and turns negative. Therefore there is a region of phase space, close to the 



point of zero recoil of the D* meson, for which the matrix elements in ( |2.4| ) are known 
directly from experiment or via heavy quark symmetry. 

It is not obvious how the form factors g{t) should be incorporated in the calculation. 
Using the matrix element given in terms of a form factor in Eq. ( p.5|) in the expression for 
the amplitude in Eq. ( p.2| ) gives non-covariant results. For a constant or linear form factor 
a covariant result is recovered if the "Z" graphs are included (ie, for a;'' < a graph with 
an intermediate state of the three particles D~, D*~^ and B^). But covariance is lost for 
more general form factors. Although the "Z" graphs are needed to recover covariance, the 
behavior of the form factor suggests that they are negligible in practice. In what follows we 
will compute in the rest frame of the decaying B~^ meson and neglect "Z" graphs. This is a 
good approximation at least when the resulting D* meson is at rest or recoiling slowly. The 
question of how to incorporate the form factor effects consistently is an open one and we 
hope to return to it in the future. 

Turning our attention to the weak matrix elements in (|2.3|) we notice immediately the 
close resemblance of these to the matrix elements in BB mixing, 

{B\Ob\B) = \flMlB^. (2.7) 

Here Ob = h^'^{l — 75)^67,^(1 — 75)^ and the right hand side parameterizes the matrix 
element in terms of Bb and the mass Mb and decay constant fs of the i?-meson. In both 
cases a heavy quark and a light quark are annihilated at a local vertex where a heavy anti- 
quark and a light quark are created. In BB mixing the heavy anti-quark is produced at 
rest. In B^ D*i~i~^ we already singled out the case of zero recoil as the one where we 
can compute the electromagnetic transition. The crudest approximation, vacuum insertion, 
works surprisingly well in the case of S-mesons, at least as indicated by calculations in 
quenched lattice QCD. Vacuum insertion gives 

{D^\OmB+) = fBfnMBMj,, (2.8) 

and the analogous result for the vector mesons. 

One may try a more sophisticated attempt at predicting this matrix element, using sym- 
metries to relate it to the measured matrix element for BB mixing. While heavy quark sym- 
metries do allow us to exchange the heavy b anti-quark for a heavy c anti-quark, isospin (or 
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SU{3), for the D*) which is needed to exchange the hght quarks does not give a relation be- 
tween the matrix elements. The operator O transforms in the reducible 3 © 1 representation 
of isospin (6 © 3 representation of SU{3)), but the corresponding operator for BB mixing 
transforms in the irreducible 3 of isospin (6 of SU{3)). We make the plausible assumption 
that the matrix element of the color octet operator Os = bT"''~f''{l — 75)mcT"7^(1 — 75)^ is 
negligible. Indeed, if one could insert a complete set of (gauge invariant) states between the 
currents in Os, gauge invariance would immediately imply that the matrix element vanishes. 
Inserting a complete set of states seems less drastic than assuming saturation by the vacuum 
alone. The assumption that Os is negligible leads to relations between the matrix elements 
of the 3 and 1 components of O. With P-j. = 1 ± 75, define 

0± = Jrf^P.ucYP^d ± Irf^P^dcYP-U (2.9) 

then and (9_ transform in the 3 and 1 representations of isospin and we have the change 
of basis relations: 



(2.10) 



The vanishing of the matrix elements of Og implies that between B~^ and D*~^ states = 
and therefore (O) = 

In order to relate the matrix elements of Ob and O we also need to make use of heavy 
quark symmetries. The first step is to write operators in a heavy quark effective theory 
(HQET) Ob and O corresponding, respectively, to the operators Ob and O in the original 
theory. If hg is an HQET field that destroys a heavy quark of four-velocity v, and hq an 
HQET field that creates a heavy antiquark with the same velocity v then 

Ob ^ Ob = 2 h{^'YP-dKl^P-d. (2.11) 

The superscript T stands for the transpose. The factor of two arises from the two ways 
of obtaining a heavy quark creation and corresponding anti-quark annihilation operators. 
Similarly 

0{x) d{x) = e"^(™''-"^=)^-^ 

X hl{x)-f^YPMx)hcixhuP-dix), (2.12) 

and analogous correspondences for Os and 0±. In the HQET there is a SU{2) x SU{2) 
flavor symmetry that rotates b ^ c and b ^ c independently. The symmetry is explicit 
between states with mass independent normalization. Using this together with suppression 
of Os and isospin symmetry we have 

{D+\d\B+) ^ 3 (D+|0+|P+) ^13 {B^IObIB') ,^ ^3. 

4 Jm^m^ 24m,' 



The factor of 1/2 in the last step is a Clebsch-Gordan coefficient. Note that the factor of 2 
appearing in Eq. ( |2.11| ) was retained in the definition of Ob so that Ob and 0+ have the 
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same normalization in the triplet tensor operator. In terms oi Bb defined in Eq. (p.7|) we 
finally obtain for states with vanishing momentum 



{D^\0{x)\B^) = e-^(— ^^flMlB^, (2.14) 

or, using f^Ms = foMi) from heavy quark symmetry, and evaluating at the origin 

(D+|O(0)|i?+) = U^M^Mj^B^. (2.15) 



Comparing with Eq. (|2.8| ) we see that the symmetry argument and the octet suppression 
assumption give the result that the B correction factor for B^ — > is precisely B^, the 
correction factor for BB mixing. 

The mixing between B*^ and D*~^ in (^.3]) is fixed by heavy quark spin symmetry M 



to have the same magnitude but opposite sign than Eq. ( p.l5|) . As a consequence of this 
sign the contributions to the rate from the two time orderings in Eq. (|2.2|) add. However, 
in the decay to a pseudoscalar, B^ D^e~^e~ , the two time orderings tend to cancel each 
other. Moreover, the electromagnetic form factors in B^ D^e^e~ are unity at t = 0, so 
the cancellation is rather good. For this reason we have not analyzed this process further. 



III. RESULTS AND DISCUSSION 

The effective Hamiltonian for the weak transition is modified by short distance QCD 
corrections, 

Kff = ^ VubV;^ {c{fi/Mw)0 + cs{fi/Mw)Os) . (3.1) 

The short distance coefficients c and cs are chosen so that their dependence on the renormal- 
ization point /i cancels the /i-dependence of operators, in such a way that matrix elements 
of the effective Hamiltonian are //-independent. Resuming the leading logs, they are given 
at /i = mfe in terms of x = {as{mb) / as{Mw))^^'^^ by c = |a;^ + ^x~^ and cs = x~^ — x^; 
numerically, with as{mb) / as{Mw) ~ 1-9 one has c ~ 1.0 and cs ~ —0.6. 

The calculation is now straightforward. The transition amplitude for B~^ D*~^e~^e~ is 
given by the vector in Eq. (|2.2|) contracted with the leptonic current and multiplied by the 
photon propagator and coupling constants. The hadronic matrix elements in Eq. ( |2.2| ) are 
computed as described in the previous section. 

There is one subtlety. For the first time ordering in Eq. ( ^.2] ) the energy denominator is 

qo + ED*+-ED = Mb -Ed, where Ed = ^ El -Ml + Ml is the energy of the intermediate 
D-meson. For the second time ordering the energy denominator is go — Eb + Eb* = —Ed* + 
Eb»i where the intermediate B* energy is Eb* = \J Efj, — M|,* + M^, . In the rest frame of 
the 5-meson and at zero D* recoil, these denominators are Mb — Md and —{Mb — Md) + 
0{1/M). For non-zero D* recoil, the denominators are different. Therefore our results do 
not depend on the form factors through the simple combination gsitB) +90(^0)) except at 
the point of zero recoil. 



6 



The differential decay rate is best written in terms of two functions that contain the 
kinematic dependence. The dependence on the vector current form factors is in 

Q{q') = Ng' i^gDitD) + Ei^E^:^rE^fB{tB)j , (3.2) 

where tx = (Mx - Mx*f - 2Mx{Ex* - Mx*) and Ng is a normahzation factor defined so 
that Q = 1 aX = q'^^^ = {Mb — Md*)^ and is approximately 

Ng^ flD + f^B^ ^I^D- (3.3) 

The second dimensionless function is a combination of amplitude and phase space depen- 
dence. It vanishes at the upper end-point as (g^ax ~ q^Y^'^- It is given by 



n,') = ^ (3.4) 

X (q' - (Ms + Moff"(q' - (Ms - Moff. 
In terms of these functions we finally obtain 

Using Eq. ( p.6|) the differential branching fraction for D*^e^e~ is, numerically, 

, I \ 2 / \ 4 

2 dT _ -I „ in-ii / rw^Vcdl \ / JbVBb \ 



tbM'— = 1.3 X 10 



dg2 ^8x10-4/ yirOMeVy 

and about a factor of 16 bigger for D*^e^e^. The integrated branching fraction, 

restricted to q^ > 0.1 GeV, is 6.3 x 10~^° assuming Q is constant, but for Q given by a single 
p-pole, Q = m2/(m2 + _ g^^^^), it is reduced to 2.2 x 10"^^ (^and, again, a factor of 16 
bigger for i?+ — > Z}*+e+e^). 



IV. CONCLUSIONS 

We have calculated the rates for B~^ —>■ D*~^e~^e~ and B~^ D*'^e~^e~. Our primary 
motivation is the observation that at the kinematic point where the D^^^ meson is recoiling 
slowly in the B rest frame the amplitude can be computed in terms of other observables, like 
the width of the D* and the B — B mixing rate. While this observation opens up exciting 
possibihties of computation and of the determination of Vub, we leave many loose ends. 
The computation makes several assumptions which need to be further tested. Among these 
the assumption that the rate is dominated by the lowest mass resonance is most suspect. 
Also, we leave unanswered a fundamental question: how can one covariantly incorporate the 
electromagnetic form factors g{t), which are here assumed to be known, into the calculation? 
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We are not discouraged by the seemingly small branching fraction obtained. The theo- 
retical questions and challenges raised (sec previous paragraph) are worth pursuing in their 
own right. Their resolution may guide us to a better method for a determination of Vub in 
exclusive decays. 

Acknowledgments This work is supported by the Department of Energy under contract 
No. DOE-FG03-97ER40546. 
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